Oles Zaburannyi 



Strong-interaction approximation for 
transfer matrix method 



June 9, 2010 



Abstract Using transfer-matrix method a correspondence between 2D classical 
spin systems {2D Ising model and six-vertex model) and ID quantum spin sys- 
tems is considered. We find the transfer matrix in two limits - in a well-known 
strong-anisotropy limit and a novel strong-interaction limit. In contrast to the usual 
strong-anisotropy approximation, within the strong-interaction approximation we 
take into account the non-commutativity of transfer-matrix components. The lat- 
ter approximation is valid for low temperatures or strong interaction in one spatial 
dimension. We observe that the Hamiltonian of the corresponding quantum chains 
contains multispin interactions. 
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1 Introduction 

Transfer matrix method allows to find thermodynamic properties for many low- 
dimensional models |[1], ||2| and shows thermodynamic equivalence for many of 
them 0, m. Important feature of the method is the possibility to establish a 
relation between thermodynamic of c/-dimensional classical models and ground 
state of [d — 1 ) -dimensional quantum models |5 1 under the assumption of strongly 
anisotropic interactions in a classical system. 

In the present paper we discuss the conditions which are imposed on a clas- 
sical system within strong-anisotropy approximation. Besides that, a new strong- 
interaction approximation which requires weaker conditions is introduced. For 
the new approximation the transfer matrix is written in the symmetric form T = 

O. Zaburannyi 

Institute for Condensed Matter Physics, National Academy of Sciences of Ukraine, 1 Svientsit- 
sldi Street, L'viv-11, 79011, Ukraine 
E-mail: zab@icmp.lviv.ua 



2 



see (|9]l, and the expansion in nested commutators similar to the Baker- 
Campbell-Hausdorff formula is performed. The quantum Hamiltonian can be found 
explicitly only if the nested commutator of an arbitrary order [X, . . . [X, 7] . . .] can 
be written in a specific form, (see lAlOl l. The strong-interaction limit reduces to the 
strong-anisotropy limit if we assume [X,Y] = 0. The strong-anisotropy limit also 
can be obtained by applying the exponential operator decomposition technique 
ISl, Q. The main difference between the operator decomposition technique and 
the strong-interaction limit is that in the latter case we must calculate nested com- 
mutators instead of assuming that oc ty. For the 2D spin-; classical Ising model 
the strong-interaction limit leads to the appearance of three-spin interactions in the 
resulting quantum chain. For the six-vertex model the new approximation gives 
the XXZ chain with four-spin interactions. 

We will use the transfer matrix method for classical 2D models fE\ with the 
total energy that can be represented as a sum over rows 

M 

E=Y^E(U^m+l). (1) 
m— 1 

where is a variable defined on the row m with IL possible values. The free 
energy of the classical 2D model (per row) can be written in the form 

where T is the transfer matrix with elements 7> e = e'^^i^t^^m+i) and the max- 
imal eigenvalue Amax which is real, unique and positive according to the Perron- 
Frobenius theorem [2] for the positive symmetric matrix T. We can consider L 
configurations as orthonormal basis {|0}^=yx °f IL-dimensional Hilbert space. 
Each configuration ^ corresponds to a base-vector in the Hilbert space. The trans- 
fer matrix T corresponds to some operator T = 1^)^^ ^'(^'l- Let us intro- 
duce the Hamiltonian of a quantum system, which is defined by the logarithm of 
the operator T 

H = -ilnT. (3) 

The ground state energy for the quantum system described by Hamiltonian OJ is 
equal to the free energy (|2]l, 

Tr I T^lnT 



Trfe-^^^H) , v / i 

''=i:^^^^7^=~-^i^^ ^^^^ =-p'^^-=f- 

Here Pq is the inverse temperature of quantum system OJ and it must be distin- 
guished from the inverse temperature of classical system ([T]| j3. We will elaborate 
approximate method which allow to find Hamiltonian ^ for 2D classical systems 
([T| at low temperatures or with strong interaction in one spatial dimension. 

The paper is organized as following. In the next (second) section the strong- 
anisotropy and the strong-interaction approximations in general for a 2D classical 
system are presented. In the third and fourth sections the methods are applied for 
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the 2D Ising model and the six-vertex model respectively. Quantum Hamiltonians 
are obtained there. Finally, in section five we compare the approximations by cal- 
culating critical temperatures for classical models. Important calculations which 
we use for the strong-interaction approximation are collected in Appendix. 



2 Approximations for transfer matrix 

In general, for an arbitrary 2D classical system with the total energy ([ij it is im- 
possible to find the transfer matrix logarithm and write down the Hamiltonian OJ 
explicitly. In this section we consider two approximate methods used to find the 
quantum system Hamiltonian Q. Both of them are based on division of energy 
E ^ into two parts, which depends on configuration ^ of single row only and 
Ey which depends on two configurations ^ and ^' of successive rows. There are 
many possible divisions, but we choose the division symmetric with respect to two 
neighboring rows 

i?(^,0 = ^+i?.(^,^') + ^. (5) 

Here ZJ, is the energy of a single row with configuration ^ and Ey[£,,£,') is the 
energy of the interaction between two rows with configurations ^ and t,' . Division 

(ig allows us to write T^^^, = T^T^ ^,T^„ where = e'P^, ^, = e-P'^yi^^^'). 
Accordingly we can rewrite T as a matrix product, T — f^fyf^^ ajjj j 
operator product, 

T= i; |^)^^^,(^'|=T-^T''T^ (6) 

where = I^^^j lO^f I is the diagonal operator and T>' = |^)r|^,(^'| 
is the off-diagonal operator. 

For further progress we have to rewrite operator T*' in a slightly different form. 
It is convenient to assume Ey{^,^) =0 (that can be always achieved by includ- 
ing nonzero value Ey{^,^) into Ex{^) or by shifting all energies by a constant 
value). That makes diagonal elements of T^' equal to unity, 7| ^ = 1- Among all 
off-diagonal elements of we distinguish the elements which are proportional 
to some parameter ty (which we will demand later to be small) and collect these 
elements in operator Y 

f if TI s, oc ty 

Y^,^'= J'^ • iif ■ (7) 
[0, m other cases . 

All another off-diagonal elements which are ^{ty) we collect in operator Y'. Sum- 
ming up 

T = e'-'-^ (1 + tyY + ^ity)Y') e'-'^, (8) 

where diagonal part is rewritten in the form = e''"'^ Until now no approx- 
imation has been made. We only assume that the parameters and ty exist. For 



4 



both approaches, i.e., strong-anisotropy limit and strong-interaction limit it will be 
later required that is small and we will neglect the terms ^(ty). For each specific 
model the individual operator Y should be constructed bearing this requirement in 
mind. 

2.1 The strong-anisotropy limit can be introduced by neglecting the terms ^{ty) 
and by a naive assumption - commutativity of X and Y in the expression for T: 

T = e'-^(l + ^T+^(f2)Y0e'-^ = e'-'V''V-'-^ + ^(f2) « e2,,x+rrY ^ 

^ ■ ^ y ' / V > / |X.Y]=0 

More precisely, the last approximation in ^ can be obtained by series expansions 
of the exponents and neglecting the terms i^(f^) and ^{ty), 

T = e'->-^ (l + fvY + ^(f2)Y') s'"^ = e^'-^+'vY + ff{t^) + ^{t^). (10) 
Finally, the Hamiltonian in the strong-anisotropy limit has the form 

H,,z = -^lnT = -^X-|Y. (11) 

For the strong-anisotropy approximation it is often assumed that tj oc ty and 
^(fy ) = 0. This conditions impose some relations on the classical system parame- 
ters ([TJ. Assumption °= ty is good to explain why approximation is called strong- 
anisotropy limit. If we consider the simplest case for model ([T]) with all in-row 
energies proportional to e^, Ex{£,) °= ex, and all inter-row interaction energies for 
different configurations not smaller than Cy, Ey{^,^') oc ey, ^ ^ then the con- 
dition tx °= ty reads 

Pe^ocQ-P'^y. (12) 

By demanding e"'^*> to be small we suppose ey to be large and e^ to be small. It 
should be noted that by demanding e"'^'"' to be small we also assume that gy > 0. 
Conditions oc ty, &{ty) = is weakest from a set of assumptions tx G ^{ty), 
^{ty ) = which allow to apply strong-anisotropy approach. 

2.2 The strong-interaction limit, in contrast to the strong-anisotropy limit, de- 
mands only ty to be small and does not put any restrictions on f^ . We can rewrite 
the expression for transfer matrix dill in the form 

T = e'-^(l +tyY + ^(f2))e'-x = c^'-.x+,y{Y+z) ^ ^^^2^ ^jg^ 

where Z is an unknown operator defined by the function Z = Z{tx\, Y). The func- 
tion Z can be expand as series in the nested commutators of operators X and Y 
(see Appendix). If a general nested commutator [X,Y]„ = [X, [X,Y]„_i] can be 
presented explicitly via some operators L, z, R (see lAlOt . all calculations can be 
performed to the very end and Z can be found explicitly. Thus, the Hamiltonian in 
the strong-interaction limit takes the form 

2t t t Q 
H„v = — -iX- -|Y- £ £ U,q■^^p.l{tx^l.c)^l,q, (14) 
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where si^pj (x) = ^ Lp=o e ^ .s^ ixe'T \ , .s/ [x) = — 1 . The strong-anisotropy 

limit (II lb follows from the strong-interaction limit (114b as expected. If we as- 
sume tx °= ty and neglect the terms &{t^) the last term in ( 1141) vanishes since 



3 Two-dimensional Ising model 

In this section both approximations will be used for the 2D classical Ising model in 
order to find the corresponding ID quantum system. The classical 2D Ising model 
is described by the Hamiltonian 

M N 

^ ^2 Jx^m.n^m.n+\ ~^ Jy^m,n^m+l^n^ (15) 

where am „ assumes two values A configuration in each row m is defined by 
a set of variables {a,„.„}^,^Y7v takes IL = 2^ possible values. We will denote 
the configurations on two successive rows by {a„} and {a^} and use the following 
notation 

N 
n=l 

N 

Ey{{an},{a:,}) = -'£jyana:,. (16) 

TheHilbert space can be spanned by the basis {|^ )}^_y^ — {^^-i |<7„)}^^ Oiv=± ' ■ 

We can establish a correspondence between the states |a), a = ±^ and the eigen- 
vectors of the spin operator s~ = Y.a \ '^)'^{'^\- Here and further on Y.a denotes 
Y,a=±^ - The diagonal part of the transfer matrix X has the form 

cTl On n=l n=l 

, N N N N 

— g j=i<yi <yNn=i «=i =e ;"=i ^ (17) 

where we omitted the direct products of the identity operators. The off-diagonal 
operator takes the form 

(^1 0]V cj[ (J^ (1=1 n=l 

In contrast to the case of diagonal operator T^, Eq. ( 117b . we can not bring the 
direct products under the exponent. To construct the operator Y, we classify all 
matrix elements T^^ } {cj',} ^^^'^ three groups depending upon a number of different 
variables in {cr,,} and {c7,^}. In the first group we collect all elements for which the 
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rows {<7„}, {<7^} are the same, i.e. the diagonal elements of T^'. Diagonal elements 

are equal to e~^, as it was discussed above, we can obtain T , , , , = 1 after 
shifting the energies by an appropriate quantity. All the matrix elements where the 
variables only on one cite differ, i.e., {a„ = ol,} ^^jn nj^j ' '^j ~ ^^^^ form the 
second group. In this case the interaction energy increases by y and the matrix 

element of T^' after energy shifting is equal to e"^ ^ . The third group consists of 
all other matrix elements for the states {g„}, {ct^} which are different on two or 

more sites. These elements are equal to ^e"'^^^ , where r is the number of sites 

with different variables a„, a^. Now, it can be seen how the small parameter fy 

should be set. If f,. = e"^ t is small, we can construct an operator Y which has all 
matrix elements equal to zero except the elements between the states which differ 
by one variable Oj (second group). All other matrix elements (third group) which 
are of order ff{ty) may be included into a non-important operator Y'. Formally it 
can be done by rewriting the sum ( 1181) as follows: 

di 0]Vn=l ,1=1 

N .7-1 N ./-I N 

+fyLL---L(8)|f^«)®l^j) (8) K)(S)i(^n\^{-Oj\ (g) ((7„I 
y=l cJi o>v«=l n=j+\ n=l n=j+\ 

I e-/^^>'({-}'{-'»(g)|a„)(g)(a;| 

n=\ n=\ 



all other 
configurations 



=l + fvY+^(f;)Y'. (19) 

(We notice that in order to calculate higher approximations with respect to op- 
erator Y' should be presented as series with respect to ty, r > 2. The term pro- 
portional to ty contains (^J) matrix elements of T^). The operator Y can be easily 
identified in terms of spin operators, 

N i-1 A' .7-1 N 

Y= LL---L(8)I^«)®I^.7> (8) k„)(g)(a„|®(-(7;| (g) (a„| 

i=\ ai On n=l n=j+l n=l n=j+l 

= 2£s}, (20) 

where = 5 Lcj (^^|- From Eqs. ( [17]) and ( l20l) we see that t^, ty, X, Y can be 
written as 

Rj J N N 

t, = ^, ty=c-Pi, X=£s)s)^i, Y=£2s^. (21) 

^ 7=1 7=1 
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In the strong-anisotropy limit we obtain the Hamihonian of the quantum Ising 
chain in transverse field 



H 



sal 



7=1 



2 -B'-y X 
— e ^ 2 « . 

j3 '• 



(22) 



To obtain the Hamiltonian in the strong-interaction limit (i.e., ty is small) we 
have to consider the commutators [X, Y]2<: 



N 

[X,Yb = 2£sj 



k>\. 



(23) 



Comparing Eqs. ( |23] | and (lAlQb we find that we have -P, ^ — 1, / — 0, Lo,i — 
i;Sy(l +4s^_iSy^j), zoj = 1, Ro,i = 1 and only one function M,o(fx) = ■s^ (tx) in 
j 

(IA12b have to be calculated. Therefore, from (114b 



H 



sil 



I As}S)^i + |e-/^^sj (l + (^) +2^ s}_is5^,) . (24) 



In comparison to (|22] | the strong-interaction approximation (l24l) implies renormal- 



ized transverse field 
of s}_iS;^s}^i type. 



1^1 



and additional three-spin interactions 



4 Six-vertex model 

In this section we will find the quantum Hamiltonian in two approximations for 
the six-vertex model \TU]. The model consists of arrows on a two-dimensional 
square lattice. Arrows are associated to each link between the nearest lattice sites. 
The total energy is the sum over the lattice vertex energies. Each vertex energy 
depends on configurations of four neighboring arrows. Moreover, only six of such 
configurations are allowed. In Fig.l the vertex configurations and the appropriate 
energies are shown. 

Not all possible arrow configurations are allowed: for each vertex the arrow 
configuration must be one of a set depicted in Fig.l by bold solid arrows. In this 
paper we will consider a lattice "rotated" by j with periodic boundary conditions 
imposed. The total energy has the form 

M N 

m=l n=l 
N 

+ ^ £'(Alm,n,Mm,n+l,Cr,„+i „,am+l_„+l), (25) 
n=l 



where the sum runs over all vertices. Arrows directions are encoded by variables 
a and /j. We will denote the variables in the two nearest rows m and m + 1 by 
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eo 









Fig. 1 All possible configurations of rows at a single vertex. The six-vertex model admits only 
the first six configurations (solid arrows) with finite energy. All configurations prohibited for six- 
vertex model (dotted arrows) have the energies that will be tend to infinity and the Boltzmann 
weights become zero. 



{a},{/j} and {a'},{/i'} (see Fig. 2). The local energy for each vertex has the 
form 

+e„5a,a2 + Sa2^^2 + 5^,' ^^2 ) + ^-^^f i^c^m ^ol + S^l S^^) , (26) 

where 5ah is the Kronecker symbol, 5^ = 1 — 5„b- All forbidden configurations 
appear in ( 126b with energies e«, which will be later sent to infinity and the Boltz- 
mann weights for the forbidden configurations will go to zero. By using the lo- 
cal energy symmetry E {1X1,112,01,02) = E{oi,02,pL\,IJ.2) and the trivial identity 

exSaia2^aiiiiSa2ti2 = ex5(j,cj2(l - ^aiMi^^" ~ ^al ^02^12 " ^al^oD we can divide 
energy in two parts 

Ex{0\,02,li\,li.2) =exda,.a2, 

Ey{o\,02,IJ.i,IJ.2) = eo5^f 5(j,^[5cj2^i2 +0 -5(1,(12 5(7,^,5(72^12 +ey5atdal da2 

+ (fioo - e^) 5(7, (72 (5(7, ^1, 5^2 + ^^1' ^(T2M2 + ^^l' ^^2 ) 

+e»5^M5^iMi5a2+5^;C)- (27) 

A reason for introducing the term proportional to zero in (127b will be seen later. 
Now, following our scheme we can write down the two components of the transfer 
matrix 

7'f(7,,}.{(7',,'} = n5(,,(,,5,„,,e-/^¥^"^>^--»+i, (28) 
«=i 



N N 

'{c7,M},{ct'.M'} 



jy ^ ^ = J^''^''e-/^£Y(f^-nfT,,+i,M«,M,,+i) J^''''™e-j3£:y(fT,;,o-,;+i,M»,M,,+i)_ (29) 



n=l n=l 



Let us consider e"^^> (°^' '^2:Mi M2) jj, niore details. As it can be seen in (|27] |. Ey de- 
pends on four variables ai , a2 , /ii , /i2 • Each of these variables can take two values 
±j, and therefore, set of four variables can take 2^* different sets of values. For 
each set of values of Oi, O2, jUi, jU2 one and only one term in (1271) is nonzero, and 
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Fig. 2 Some permitted configurations for the six-vertex model. We show two neighbor rows to 
explain a definition of variables a and ju. 



therefore we can write 

+ e-^-5^f + 8^l8a,^,). (30) 

We are free to redefine the energies, e^ — cq, ey Cy — eo, eo — > 0, and to 
send the energies of forbidden configurations to infinity, e"^''°° — >■ 0. By denoting 

ty — e'l^^y we get 

^-PMaua„^^u^^2) = 5,,,^^ 5^^^^ +tySSf8^l8^^. (31) 

Now we have 

N 



n=l 



n=l 

N N N n+i 

n=l j=ln=l 



5;^-5j52;;5<^^5<^,,,,,.,ifjeven ^^^^^ 



5^^+ ' 5^' 5^'^ S„,, .d„, „ . , , if y odd 



Using the identities Sab^bc = SatSac, SabS§ = SabSa we can rewrite in the form 

N 



n=\ 



;=in=i 



(33) 
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The matrices T, and T'' are of the size 2'^" x 2"^™ and the quantum Hamihonian 
H = — 4 In T should be defined in the 2^^^ -dimensional Hilbert space. We can show 



that in the particular case ( l29l) the space dimension can be reduced to 2'^. Let us 
recall that we are looking for thermodynamic properties of the system described 
by the partition function TtT^ (|2]l. Matrix ( |29l l can be represented in the form 
j,y _ J,} _^ + ^(f^) and by using the trace properties we may write 



TrT 



M 



-tyTl 



M 



= Tr 



((r 



x\2rpy^ 



f ^(f2))r 



(34) 



The trace over space {cr,/!} can be splitted into two partial traces over {a} and 
over {pi}. From Eqs. (128b . ( |29l l one can easy write down the matrix elements for 

[{T'^ fTl)^ and Mty{{T'fTy 
result is as follows 



rx\2Ty\^ 1 j^y ^jjj perform the partial trace Tr{^}. The 



Trr^ =Tr{^}Tr{^} 



{{r)^Tlf +Mty{{T 



=Tr{^}e"''''^'^^'<'^'^.+' 



N N N n^j 

n=\ j=\n=\ 



+ ^(f2). 



(35) 



The latter expression in Eq. (1351) is Tr{fj}(f'^f''T^)''^ + ^{ty), where the matrices 
T'', T^' are defined by their elements 



n=l 
N 



(36) 



We can look for the ground-state energy of Hamiltonian H, which acts in the 
2^^ -dimensional Hilbert space {a} and is defined by the logarithm of the matrix 
T = f -*7>'f The reason that we can perform space dimension reduction is as 
follows: For any two states {cr}, {cr'}, which differ by only one Oj, there is only 
one configuration {/j}, for which the matrix element of T remains non-vanishing 
in the limit &{ty) = 0. From here we will omit the tilde and use the notations T, 
and T^' to denote the matrices that depend only on {a} variables. Matrices 
T^, ( 1361) are represented in the form suitable for obtaining operators X and Y 
mentioned in ([8]l. 

N N N 

x = II(8)l^«)V.+.(8)(^«l' 

7=l{cj}«=l n=l 



N N N n+j , „/ A' 

y=l {CT,(7'}«=1 n=\ ' n=\ 



(37) 
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If the states defined by {a} are considered as the spin states |^) = |-^) = 
we can immediately recognize the spin operators 

L |aj) ® \aj+i)5a.a.^, {aj\ {oj+i \ = 2(8)8)^^ + 1/4) 

(38) 

I I |cTy)®|ay+i)5^5^+'5;/-'(a;|®(a;+i|=2(s5s5+i+s>^^i), 

(39) 

We have almost all that we need to construct the quantum Hamiltonian in both 
approximations. At last we have only to recall Eq. ( 128b to define tx and collect 
essential variables and operators 

h = -^. ty = ^-^'^\ X = 2£s^s)+i + i/4, Y = 2£s>^+i+s>;;^i. (40) 

^ ./=i ■ ■ y=i 

In the strong-anisotropy limit dl lb we can immediately write the quantum Hamil- 
tonian, which corresponds to the spin- 1/2 XXZ chain 

i=i 

where we omit the insignificant constants. For the strong-interaction limit we have 
to calculate the commutator 

[X,Y]2, = 2^^ £ (1 -4s5_iS}+2) (sjs;j+i (42) 

This commutator has the form dAlOb withP, g = 1, / = 0, Lo,i = ~4s^_iS^^2) 

•(s^'Sy^j +s-ySv|^j), Zo.i = 2, Ro,i = 1. Since P = 1, we need calculate only one 

function in ( IA12b si^Q,i{tx) = ^{tx) = ■s^{—tx), in order to construct operator Z 
( IA13b . The Hamiltonian ( 1141) becomes 

H,, = 2 £ .,s)s}^i - yf^^y (1 + (1 -4s5_,s}^2)) (s>5+, +s^s^+,). 

(43) 



5 Conclusions 

To compare the results provided by both approximations we will discuss the criti- 
cal temperature for the Ising model which corresponds to quantum phase transition 
of the Hamiltonians obtained in the strong-anisotropy limit (122b and in the strong- 
interaction limit (1241) . Both quantum (1221) and (124b models are particular cases of 
generalized spin-^ XY chain for which a critical point can be found by Jordan- 
Wigner and Bogolubov transformations [lU . Moreover, both quantum models 



12 



and ( |24] | have the critical point at the parameters which correspond to the 
same value for the critical temperature of the classical 2D Ising model 

PsalJx PsalJy too r a a\ 

exp = 1 , j3.„7 = j3,oz • (44) 

The fact that more accurate approximation, i.e. the strong-interaction approx- 
imation, does not give any improvement for the value of the critical temperature 
can be explained based on analysis of the well-known exact equation for the Ising 
model critical temperature, 

smh smh = 1 . (45) 

2 2 

Eq. (l45T l can be rewritten in terms of t^, f,, (1211 1 as follows: 

t, = arsinh (^YTf? ) = + ^^^'^ ' ^"^^^ 

In the strong-interaction limit i^{ty) = and we have tx °= f,, that is the requirement 
of the strong-anisotropy limit. The fact that both gave the same critical temperature 
can be explained in terms of interactions J^, Jy and inverse temperature j3. In space 
of parameters Jx, Jy, P the strong-interaction limit covers much wider region than 
the strong-anisotropy limit. But near to the surface of critical temperatures defined 
by Eq. ( 145b these regions coincide. More precisely intersections of two surfaces 

given by exact solutions to ( 1451) and by the strong-interaction limit e"^^ =0 is 

a line defined by the strong-anisotropy limit e'P^ =0, ^ = 0. This transparent 
geometrical interpretation unfortunately can not be simply depicted because all 
regions which are object of our interest are infinitely distant. 

For the six-vertex model the critical temperature in the strong-anisotropy limit 
can be calculated from isotropy condition for the XXZ Hamiltonian, (141b 

o-P'^'y = Pexex- (47) 

The critical temperature in the strong-interaction limit requires a study of the 
quantum spin chain described by Hamiltonian ( 1431) for which no exact results are 
available. From the exact equation for the six-vertex model critical temperature, 

e-fe^-v^e-fe^-, = 1^ (48) 

we can draw conclusions similar to those derived for the Ising model. In fact, 
equation ( 148b can be rewritten (see (140b ) as 

tx=^^^^ = ~j + ^{t^), (49) 

wherefrom the equivalence of two approximations near the critical temperature 
may be expected. In Fig. 3 the results for critical temperature given by strong- 
anisotropy approximation and strong-interaction approximation are shown. We 
can see that the strong-anisotropy limit gives reasonable results even when the 
classical system become isotropic. 
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1 3 a 



Fig. 3 Dependence of relative inverse critical temperature ji on spatial anisotropy a. For the 
Ising model ^ = fij^, a = y^, solid line coiTesponds to exact solution <45t . dashed line cor- 
responds to the strong-anisotropy and strong-interaction limits l l44t . For the six-vertex model 
P = Pcx, a = dash-doted line corresponds to exact solution l l48l l, doted line coixesponds to 
the strong-anisotropy limit {47). 



In spite the fact that the novel strong-interaction approximation does not im- 
prove the critical temperature of the considered classical 2D models, the strong- 
interaction approximation have several advantages. First, we can study isotropic 
classical systems. Second, we can apply strong-interaction approximation as low- 
temperature approximation for classical systems with arbitrary interactions. In- 
deed we can achieve ty °^ e'l^'^y to be small in two ways: by assuming strong inter- 
actions along one direction e,. — > oo or putting low temperature 15^°°. 



Appendix: Representation of exponent products via nested conunutators 

We will start from the expression for transfer matrix <13K where an unknown operator Z was 
introduced e''^(l +tyY + iff{t^))e'''^ = e2''^+'>(^+^) -I- ^(f^). By expanding the right-hand side 
of the latter equation in powers of ty up to first order we derive for the linear over ty terms the 
following result 

e',XYe',x= /' rfTe2°'X(Y-hZ)e-2"-.Xe2'.x. (a1) 
Jo 

The left and the right multiplication by e^'''' and the substitution t — > /q dr — > j j\dT: 
lead us to equation: 

„v,.. _„v. , /■' ;^ T«f«(fx,Y]„-h[x,z]„)^ 

dT, 



Y=i/" e"'-X(Y + Z)e-"-XdT=i / £^ 

-'-1 ■'-^11=0 



[A,B]„ = [A,[A,B]„_|], [A,B]o = B. (A2) 

Now we can integrate over t 



Y = 3 I ^."([X,Y]>. + [X,Z]„) = £ tl 



(2n+\)\ 



(A3) 
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Let us presume that operator Z can be expanded in terms of nested commutators with unknown 
coefficients A* : Z = £ ^ ?* [X, Y] ^ . We can show that the function Z(f ^ X, Y) is even with respect 
to the first argument, 



grrXg^Ygd-X-j-l^ |^g2(,X+/j(Y+Z(/,X,Y)) \ ^g-2/,X-rv(Y+Z(r,X.Y)) 



(-g/jXg/.Yg/jX-j - 1 ^ g-fvXg-r,,Yg-/,X ^ g-2/,X-r,,(Y+Z(-r,-X.Y)) 



Z(«,-X,Y) =Z(-f,-X,Y), 
(A4) 



that means that all odd A2k+i are zero Aik+i = 0. Thus we will look for the following series for 
Z: 

7^=i^t?[X,YU (A5) 

k=i 

Substituting l lASt into iA3\ and reordering the sum (Cauchv product) £ L^("i^)=L T, 

n=Ok=0 1=0, /=0 

7) we can write 

Y = E , ^ , f?'[X.Y]2,- + f (y 1 /?'[X, Y]2,-. (A6) 

Equating the coefficients at equal powers of f, we have the expressions for A211 



en(2/-27)!(2/+l)!^f2/+l)! ' ^ ' 



;t'o(2'--2i)!(2j+l)! (2/+1)! 
From this equations we can write a recursive representation for Aj., 

where f?'.^',) = ■ -i, . It should be mentioned that for the Bernoulli numbers similar 

\2j+lJ (2;+l)!(2i-2;)! 

recursion representation exists, 



B2(i-j) + 5, ! > 1 Bo = = -5- (A9) 



Our new task is to find the expansion l lASt . Evidently it is impossible to do for arbitrary op- 
erators X and Y. We will consider only the case in which the commutator [X, Y\2k have following 
periodical (with period P) structure: 

Q 

[X, Y]2/>,+2/ = £ U.q (z,,,) R,,„ (AlO) 

q=l 



where P,2,r 6 N, / = 0,P— 1 and the operators L;.,^ = L/^p.^y, R; = R;+p,^, z/ ,j = z;+p,j. For- 
mula <A10t links nested commutators of the orders 2Pr + 2/, r = l,co by 2 rules. For each 
/ = 0,P — 1 this rules can be different. With increasing the commutator order by 2P the expres- 
sion for commutator is multiplied by operators z with cons tant s andwich multiplication by the 
operators L, R. If commutator [X,Y]2Pr+2; have the form l lAlOt starting from / > 1, we can 
separate in l|A5) the terms k= 1 , r' — 1 and perform all computations for the redefined operator 
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Z. Further progress is possible due to fact tiiat an exponential generating function si {£) for the 
coefficients A2* can be suggested 

^ ' 2k\ sinhx ^ ' 

From Eq. l lAl It we can find a generating function for s/pj (x), 

^^,(;,) y ^!ll±^^Pr+2i ^ 1 y e-^^ U^) (A12) 

Finally, the unknown operator Z takes the form 

Q P-l 

Z = £ E L,,,^p,,(f,z,,,)R,.,. (A13) 

9=1 /=o 

Having obtained this operator we can take logarithm of the transfer matrix <13t and obtain 
quantum Hamiltonian in the strong-interaction limit l fT4T l 
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